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The scaling properties of the roughness of surfaces grown by two different processes randomly 
alternating in time, are addressed. The duration of each application of the two primary processes 
is assumed to be independently drawn from given distribution functions. We analytically address 
processes in which the two primary processes are linear and extend the conclusions to non-linear 
processes as well. The growth scaling exponent of the average roughness with the number of appli- 
cations is found to be determined by the long time tail of the distribution functions. For processes in 
which both mean application times are finite, the scaling behavior follows that of the corresponding 
cyclical process in which the uniform application time of each primary process is given by its mean. 
If the distribution functions decay with a small enough power law for the mean application times 
to diverge, the growth exponent is found to depend continuously on this power law exponent. In 
contrast, the roughness exponent does not depend on the timing of the applications. The analytical 
results are supported by numerical simulations of various pairs of primary processes and with dif- 
ferent distribution functions. Self-affine surfaces grown by two randomly alternating processes are 
common in nature (e.g., due to randomly changing weather conditions) and in man-made devices 
such as rechargeable batteries. 

I. INTRODUCTION 

A. Motivation and background 

The scaling approach to kinetic roughening of surfaces has been found to be very useful in a large number of 
systems [1-7]. It is based on their self-affine property, reflected by the power-law dependence of the roughness on the 
system size in the, steady state, large-time regime. In the early regime, the roughness also increases as a power law 
of the time. Systems in which such a behavior of the surface roughness was observed and analyzed include different 
(vapor, electrochemical, epitaxial, etc.) deposition processes, burning front, malignant tumors, and others. Similarly, 
self-affine rough surfaces are formed in opposite processes of recession, desorption, dissolution, etc. 

Very recently we have addressed the scaling behavior of cyclical processes in which two different processes (typically 
one growth and one recession) follow each other regularly. Such cyclical processes are abundant in natural and 
artificial systems. Many natural growth processes are affected by changing degree of illuminations through the 
daily cycle or seasonal variations in weather conditions. Rechargeable batteries, charged and discharged periodically, 
provide an example of a practical application in which a rough metal surface accumulates on one of the electrode 
(and cause a short when it reaches the other electrode) . Chcmotherapeutic (or radiation) treatment of a malignant 
tumor applied periodically and the expansion/curtailment of a bacterial colony depending on the availability of 
the nutrient, may provide examples in the biosciences. Experiments on the roughness of silver surface grown by 
cyclical electrodeposition/dissolution were found to be consistent with the scaling predictions. The main results of 
our theoretical investigations are summarized in Sec. II. 

While studying these cyclical processes it was also realized that many of them are not exactly periodic. Namely 
the duration of one (or both) primary process, may not be uniform. For example, the charge and discharge time 
of a battery maybe very irregular depending on the way it is utilized. If a growth process depends on atmospheric 
conditions such as alternating rain and sunshine, the length of the application time of each primary process will 
also fluctuate drastically. That may be the case for alternating dry/wet erosion processes of a natural rock bed or a 
5^ , man-built wall. 

So two crucial questions which ought to be addressed are: To what extent the cyclical theory applies to systems 
with non-uniform application times of the primary processes? And if the cyclical theory fails to describe them, how 
are their scaling properties modified? The goal of the present work is to seek and provide the answers to these two 
essential questions. 
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B. Basic scaling concepts and relations 



The roughness of sclf-affine surfaces created by growth or recession processes have been described using scaling 
concepts [8]. The width W(L,t) of a surface linear dimension L at a time t is given by 

W(L,t) = ((h(r,t)-(h(r,t))) 2 y, (1) 

where the angular bra-kct (} denotes average over both lateral sites and different realizations of surface configura- 
tions. W(L,t) was found to scale as [8] 

W(L,t)~L a g(L/£(t)), (2) 

where ~ i 1 / 2 is the lateral correlation length, g is the scaling function which behaves such that for large time 
(t 3> L z ) W ~ L a (a is the roughness exponent), while in the early-time regime (t <C L z ), W <~ t 13 (where (3 = a/z 
is the growth exponent). Different surfaces are classified into universality classes which share the same set of scaling 
exponents. Some generic universality classes are mentioned in subsection C. 

For rough surfaces formed by cyclical processes a suitable scaling law in terms of the number of cycles n was 
introduced [9] [10]. The width W(L,n) for a system of linear size L obeys the following scaling form: 

W(L,t) = L a g c (L/Un)), (3) 

where £ c (n) ~ v}l z is the lateral correlation length and g c is the cyclical scaling function. 
For large time b> L z : W ~ L a , 

while for n <C L z : W ~ rfi , 

where = a/z is the growth exponent. Again, different universality classes can be defined depending on the values of 
the exponents a and 0. In most of the cases studied by us, the scaling exponents of the cyclical process were identical 
to those of one of the primary processes (the so-called dominating process). 

In our previous papers on cyclical growth [9] [10] the durations of the primary processes (in every cycle) were 
assumed to be uniform. However, in many realistic situations the durations of the primary processes will vary from 
cycle to cycle. 



C. Growth Models and Universality Classes 



The growth processes fall into different universality classes [2] . All processes within one class share the same expo- 
nents and their asymptotic continuum stochastic equations differ at most by irrelevant terms (in the renormalization 
group (RG) sense). Using the symbolic index i = 1,2, .. to denote different processes, the ones we consider here follow 
growth equations of the form: 

^A=A i {h} + r li {r,t)+v i , (4) 

where Ai{h} is a local functional depending on the spatial derivatives of h(r, t) and the noise f]i(f, t) reflects the the 
random fluctuations in the deposition process and satisfies 



(Vi(r,t))=0, (5) 

and 

(Vi{ri,ti)vi(r2,t 2 )) = 2D i 6 3 {f 1 -f 2 )6{t 1 -t 2 ), (6) 

where d — d — 1 is the substrate dimension. 

Some of the generic growth processes we will consider are: 

(i) Linear: Random deposition (RD), Edwards- Wilkinson (EW) [11] and Mullins-Herring (MH) [12] [13] [14]. 

(ii) Nonlinear: Kardar-Parisi-Zhang (KPZ) [15] and Molecular Beam Epitaxy (MBE) [16] [17]. 
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D. The structure factor (SF) 



The structure factor is defined as S(q,t) = (h(q,t) h(—q,t)), where h(q,t) is the Fourier transform of the height 
h(f,t). In the theoretical analysis of surface growth it is convenient to work in Fourier space and compute first the 
SF rather than the roughness W itself. Also, in some experiments the surface is probed by scattering processes which 
provide its structure factor. 

Dynamic Scaling hypothesis Eq. (2) can be translated to the Fourier space such that [2] 

S(q,t) = q- g - 2a g{q/t- l / z ). (7) 

Surface width W can be readily calculated from S(q,t) using the relation W 2 (L,t) = (l/L d )^2~S(q,t). 

The organization of the rest of the paper is as follows: In the next section (II) we briefly review the theoretical 
analysis of cyclical processes and their main results. In Section III Some exact results are derived for the case of two 
linear primary processes with non-uniform application times. In Section IV various probability distributions of the 
random application times of the primary processes are considered. Results of numerical simulations for both linear 
and nonlinear primary processes are presented. Sec. VI contains a summary of our results. 



II. REVIEW OF CYCLICAL GROWTH PROCESSES 



Any cyclical growth process consist of two primary processes. The durations for the first and the second processes 
are T\ = pT and T 2 = (1 — p)T, respectively. The period of one cycle is T = T\ + T 2 . The cyclic growth equation 
(equivalent to Eq. (4) of a simple growth/recession) in terms of the basic two processes is 

dh 

— = [ ai h + m + v 2 ]e( P - f(tj) + [a 2 h + m + v 2 ]e(f(t) - P ), (8) 

where f(t) is defined as the fractional part oft/T and Q(x) is the unit step function. 



A. Linear Primary Processes 

When both the primary processes are linear and their durations (Ti and Ti) are constants the above Eq. (8) can be 
solved exactly to yield the SF. These results are is briefly described below and details may be found in our previous 
papers [9] [10]. 

For a linear process Langevin equation of the form of Eq. (4) can be easily solved in Fourier space to yield the SF 



S(q,t) = cxp{-2a(q)t}S(q,0) + ^[1 - exp(-2a{q)t)}, 



(9) 



where S(q, 0) is the SF at t = which contains the information of the initial roughness. We assume S(q,0) to be 
zero, i.e. the growth starts from a flat substrate. This will not affect the asymptotic scaling exponents. 

We index the two primary processes by i = 1,2 respectively and define ai = ai{q)Ti. During the 1st cycle of the 
cyclical growth, the structure factor generated by the first primary process (of duration T\ = pT) is assigned as the 
initial condition for the second primary process. The second process lasts for T 2 = (1 — p)T to yield the structure 
factor S c (q, 1) of the cyclical process after the first cycle. This is again used as the initial structure factor for the first 
process in the 2nd cycle. If we keep iterating like this, the SF after n cycles becomes 



S c {q,n) = 



— exp(-2a 2 T 2 ) (1 - exp(-2oiTi)) 

Oi 

+ — (1 - exp(-2a 2 T 2 )) 

«2 



1 — exp(— 2a c n) 
1 — exp(— 2a c ) 



(10) 



where a c = a c T with 



a c = [a 1 p + a 2 (l - p)}- 
In the scaling limit of small q, Eq. () for the cyclical SF reduces to 



(11) 
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S c {q, n) [1 - exp ( - 2a c (q) T n)] , 



(12) 



where the effective noise strength for the cyclic process is defined as 

D c = pD 1 + (1 - p)D 2 . (13) 

In terms of the effective cyclic parameters a c and D c , the SF (12) for the cyclical process looks very similar to that 
of a generic linear growth process (see Eq. (9)) if the time variable is replaced by the number of cycles n. Hence, the 
dynamic scaling analysis (Eq. (7)) of the structure factor can be used to determine the scaling exponents in the case 
of a cyclical growth. 

In the regime of saturated roughness, after large number of cycles n — > oo, Eq. (12) yields S c (q,n) ~ a D ^ . The 
roughness exponent of the cyclic process is determined by the q — > divergence of . Since a c — [aip + a 2 (l — p)]T 
and di(q) ~ q Zi , it is the process with smaller z% which dominates the asymptotic cyclical roughness. In the growing 
phase of the interface roughness, the number of cycles n is multiplied by a c (q) in Eq. (12). The process with smaller 
Zi will again dominate in the q — > asymptotic limit. Therefore, the primary process with the smaller dynamic 
exponent carries over its scaling exponents to the combined cyclical process. 



B. Non-linear primary processes 



We extended our analysis to non-linear processes as well by using an approximate RG approach. In this approach 
we "set aside" the non-linear terms for the first few RG iterations in which all fluctuations on time scales smaller 
than the time one full cycle are integrated out. These few initial steps deals only with the linear parts of the two 
primary processes and result in an effective linear process, as described above. At this stage the non-linear terms 
are added back to the stochastic equation with their bare couplings multiplied by p (respectively, 1 — p) to take into 
account the fraction of time of their operation. The next iterations of the RG process then may be applied in the 
usual manner. This approach assumes that the approximation made in the few initial iterations will not change the 
ultimate fixed-points to which the RG flow take the system. We expect this to be the case if the system location in 
parameter space is not too close to a separatrix between two basins of attraction. 



III. ANALYTICAL CONSIDERATIONS OF RANDOMLY ALTERNATING PROCESSES 



A. Linear primary processes 



Now we are ready to consider the more general case of nonuniform durations of the two primary processes. The 

• • (k) (k) 

time elapsed during the first and the second processes in the k application are denoted by t\ and t 2 respectively. 
So t[ k) and 4 fe) 

are assumed to be random variables obeying specific probability distribution functions (pdf) . We also 

define two new random variables = exp[— 2a\{q)t^] and rjk — exp[— 2a 2 (q)4^L in terms of which the structure 
factor of the cyclic process will be expressed. 

The SF, after the first primary process has taken place in the first cycle, turns out to be 

S(qM iy ) = -^t[1 - exp(-2a 1 (g)t)] = -^(1 - (14) 
ai(q) ai{q) 

where S(q, 0) — is assumed for simplicity. Hence, after one complete cycle the SF becomes 

S c (q, 1) = (^)m(l - £i) + (^)(1 - t?!), (15) 
ai a,2 

which can now be used as the initial SF for the first primary process of the next cycle. After the end of the 2 nd 
cycle, the SF is 

S c (q, 2) = (^i) {&»7iffc(l - + %(1 - £ 2 )} + (^) u 2 %(i - m) + (i - m)} ■ (16) 

Cl\ 02 

Proceeding in this manner we finally arrive at the general expression of the SF after n cycles 
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j-*. I n—1 n n 

S c (q,n) = (^) II ^+iU^ 1 -^)+Vn(l-in) 



r—1 i— r+1 j—r 



n—1 n 



+ (^) \ E( II & II - Vr) + (1 - r?n) \ (17) 

(Note the similarity of this sum of products of random variables to that encountered in ID diffusion in random 
environment [18]. For that system a similar expression yields the first passage time and the random terms in the 
products are the ratios of the transition probabilities in the two directions between neighboring sites). Suppose the 

duration of the primary processes and are distributed with the pdf P{t\) and Pfa), respectively. We need 
to take the average with respect to these P{ti) of the SF in order to obtain the scaling behavior of the roughness. 
However, the above expression of the SF contains various products of the random variables £ and rj and in such 
situations the average behavior might differ from the typical behavior. Since by definition = exp[— 2ai(q)ti] and 
i]k = exp[— 2a 2 (g)t|], they always take values less than one. Thus, (£ m ) < 1 and {i] m ) < 1 (m = 1,2....), i.e. all the 
moments including the averages remain less than unity, and also In (e~ 2ait i ) < 0. This is valid for any normalizable 
probability distribution due to the fact that the durations are always positive (t* > 0). As a result, the SF is not 
dominated by the higher products of and r]k , rather the most dominant contribution in the sum will come from the 
terms linear in the random variables. In that case, the fluctuation around the average will be small and the typical 
behavior will be well represented by the average. By taking the average of the Eq. (17) 

(S c (q,n)) = (^)/>P I 11 ] (l-&)+Vn(l-Zn) 

* r=l \ i=r+l j=r 

+ (^)(E II & II Vj)(l-Vr) + (1-Vn) 
\ r—1 V i— r+1 j—r-\-l 



. , I n — 1 n n 

= (-^-) E( n ® n ^))(! - + (»?xi - (0) 

I r—1 i—r-\-l j—r 



n—1 n 



+(^) \ E( II <0 II MX 1 - <»?» + (! - !> > ( 18 ) 

I r=l i=r+l j=r+l \ 

where (£) = (cxp[— 2ait\]) and (ry) = (cxp[— 2a 2 i 2 ])- The above expression (18) of the SF has the form of a geometric 
series which can be easily summed with the result 



■A (<!■»)) = \( — )(v)(l (0) + ( — )(1 - (V)) 



(19) 



This form of the average SF resembles the SF ( Eq. (10)) for the case of uniform durations of the primary processes. 
In fact, one can easily arrive at Eq. (10) starting from Eq.(19), once the uniformity of t\ and t\ are assumed, i.e. 
P(t) = S(t-T). 

The scaling exponents in this case of random application times may be extracted using the small q divergence (as 
done in the uniform case also) of the average SF (19), given in terms of (£) and (77). The averages of the variables £ 
and 77 have been defined as 

(0 = J e- 2a ^P(h) 

= A(« - 2ai), (20) 



and 



(V) = J e- 2a ^P(t 2 ), 

= P 2 (s = 2a 2 ), (21) 
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where -Pi(s) is the generating function (Laplace transform) of Pi(t). Laplace transforms of similar random time 
intervals were considered by Godreche and Luck [19] in the context of renewal processes. In such stochastic processes 
events occur at the random epochs of time ti,t 2 --- from some time origin t — 0. Change of sign of the position of 
a random walker is one such example. The intervals of time between those events are treated as independent and 
identically distributed random variables, very similar to the random duration of the primary processes in our case. 
Depending on whether the mean of the probability distribution P(t) is finite or not, two different types of behavior 
are expected, as already observed in continuous fractal-time random walk [20] and in the case of renewal processes 
[19]. 

(i) Random application times with finite means: 

When the mean of the probability distribution is finite, we can write the averages of the random variables £ and r\ 
in terms of a series 

(0 = 1- 2ai(fi) + higher terms, (22) 

and 

(77) = 1 - 2a 2 (t 2 ) + higher terms. (23) 

In the scaling limit of q — > 0, i.e. small a i} the higher order terms can be neglected with the result 

(0 = l-v^ih) and (77) = l-v 2 q Z2 {t 2 ). (24) 

Again invoking the asymptotic limit q — ► 0, the primary process with smaller z% (dynamic exponent) dominates 
the average SF and hence the scaling behavior of the roughness. Therefore, as long as the probability distribution 
has a finite mean, the scaling exponents corresponding to the average roughness are identical to those of the uniform 
duration case, with the average duration being the uniform period. In the next section, this will be shown explicitly 
for various probability distributions. We will also present results of numerical simulations which shows the validity of 
our conclusion even if one or both of the primary processes contain nonlinearity. 

(ii) Random application times with diverging mean: 

In some specific cases of broad distributions the average of the probability distribution diverges. One such example 
will be distribution with a fat power law tail 7j~(^ +1 ) with /1 < 1 (this is similar to the behavior found in continuous- 
time random walks, giving rise to the so-called fractal time behavior [20]. In this case, (£) and (77) in the scaling limit 
(q — > 0) are given by [19] 

(f)^l-/si(2ai)" and (r?) ~ 1 - K 2 (2a 2 )". (25) 

By putting a, ~ q Zi , we obtain 

(0 ~ 1 - /site)"* 1 and (77) ~ l-«£(g)"* 2 . (26) 

Hence, the dynamic exponents are modified as fiZi. The process with smaller z% will still dominate and yield the 
cyclical dynamic exponent z new = fi min{z\, z 2 }. Physically this modification is the consequence of the modified 
relation between the time and the number of cycles. If we ask what is the total time elapsed after n applications of 
the two primary processes the answer is as follows: As long as the mean duration of one application is finite (i.e., 
H > 1 for distributions with an algebraic tail), the total time is proportional to n times this mean duration. For 
power-law distributions with /i < 1 and diverging mean, the dominating contribution to the total time comes from a 
finite number of applications, each of them lasting a time of the order of n 1 ^. This also implies that for such systems 
in which the roughness grows with time as t 13 (where is the growth exponent of the primary process with the smaller 
z), it will grow as n 13 ^ with the number of applications n. 

After a large number of cycles, once the roughness becomes saturated, the roughness exponent a of the grown 
sclf-affinc surface should not depend on the random application times of the primary processes. As a result, the 
exponent a will remain unchanged from the corresponding uniform duration case, which is nothing but the a of the 
primary process with smaller z. This is consistent with our prior observation that the growth exponent is modified 
as Pnew — 0/t 1 , so that the scaling relation ^ - = z new remains valid. So for /j < 1, the exponent changes 
continuously depending on the values of ji. 

For two primary processes with durations having two different values of fi of the power law distributions the new 
dynamic exponent becomes z new = min{^\Z\, /J, 2 z 2 }, and = a/z new , where a — j{z new — (d — 1)}. If one of the 
processes (say for i = 1) has a finite mean while that of the second (i = 2) is diverging, the value of the z new will be: 
Znew = min{zi, fj, 2 z 2 }, with similar expressions for a and 0. 
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B. Non-linear primary processes 



(i) Application times with finite mean: 

In this case we expect our approximate RG approach, introduced for periodic processes and described in the 
previous section, to still be applicable. The few initial RG iterations will coarse grain all time-fluctuations smaller 
than a typical time of the order of the average time period (the sum of the mean application times of the two primary 
processes). Although the distributions will always allow remaining fluctuations from cycles longer than the mean, 
their probability decreases with their length and we do not anticipate them to be relevant to the long time behavior. 

(ii) Random application times with diverging mean: 

In this case the approximate RG cannot be applied since the dominating effect is contributed from a few extremly 
long applications, rather than from the accumulation of n applications. This behavior strongly suggests that each of 
the primary processes should be renormalizcd independently form the other and their relative contributions should be 
compared based on their fully-renormalized scaling behavior (contrary to the case (i) above where the renormalization 
is carried on after the two processes have been combined to a single effective process). 

If that would be the case, the former result z new = min{iiiZi, ^2^2}, will hold in general with the z(i) of the two 
primary processes whether they are linear or not. The exponent /3 new will be given by the larger of (3{i)/n(i). a will 
simply be that of the dominating process determining z new . 

(iii) Random application times with one finite and one diverging mean: 

If one of the two processes (say the first one) has a finite mean application time while it diverges for the other (the 
2 nd ), some special care is required. We can again imagine performing a few initial RG iterations until the average 
duration of the first process is of order one. The duration of the second process will be rescaled by the same rescaling 
factor but will continue to be dominated by the few extremly long applications. If we have z-i < z\ to begin with, 
then z new = [i2 z 2 as before. If, on the other hand we begin with Z\ < z%, then z\ is unaffected during the initial 
RG iterations and we have to compare the contribution from n applications of the first process with the few extremly 
long applications of the second one. Again, we expect the smaller of z\ and ^2^2 will dominate and yields z new of the 
combined process. In short these results amount to assigning the value [i\ = 1 for the finite-mean process and using 
again the same expressions for z ne w and /3 new in (ii) above. 

IV. NUMERICAL SIMULATIONS OF RANDOMLY ALTERNATING GROWTH WITH DIFFERENT 

APPLICATION TIME DISTRIBUTIONS 

Below we will consider different probability distributions P(t) (both discrete and continuous), and try to determine 
the scaling behavior using the SF (19). The variables defined in the previous section will continue to be used. The 
average duration of the two primary processes will be denoted by Ti and T 2 for all types of probability distributions 
(as long as the average exists and finite). 

We have carried out numerical simulation of various linear and nonlinear discrete growth models to verify our results 
corresponding to different probability distributions. The system size used in our simulation was varied between 128 
to 4096 lattice units. Periodic boundary condition is employed to keep the finite size effects to a minimum. A typical 
cycle consisted of average deposition (Ti) of 8 layers (average number of particles per site) and the same amount of 
average desorption (T2). The maximum number of cycles n was changed between 512-8192 to reach the saturated 
roughness phase. The roughness data was taken for ~ 1000 — 6000 independent runs to average over the realizations 
of the random deposition noise as well as the probability distributions for the random duration. 



The durations of the primary processes are equally probable between two limits T max and T m i n . Hence, the 
probability density is P(t) = ^ — d ^ T . This yields 



A. Uniform Distribution 



(0 = (exp(-2ai*i)) 



(27) 




(28) 



1 



[exp(-2aiT r 



) 



- exp(-2aiT mo J 



(29) 



2ai (T maa: 



min 
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exp(-aiT max ) exp(-aiT min ) 



2ai(T„ 



Train) 



exp(oiT maa; ) exp(-aiT mm ) - exp(-aiT maa; ) exp(aiT mm ) 



(30) 



= cxp(— 2a\T\ 



sinh{2a 1 (AT 1 /2)} 



{2a 1 (AT 1 /2)} ' 

where 7\ = (T min + T max )/2 is the mean time and ATi = T„ laa; - T min . Similarly, 

sinh{2a 2 (AT 2 /2)} 



(■q) = cxp(-2a 2 T 2 )- 



{2a 2 (AT 2 /2)} 



Clearly, the uniform duration case is recovered when ATi vanishes. 
In the scaling limit of small q, 



since 



sinh{2a.i(ATi/2)} 
{2a,(AT,/2)} 



(£} <~ exp(-2aiTi) and < r? >~ cxp(-2a 2 T 2 ) 
1 as a, = g Zi — > 0. In this limit, the SF reduces to 



S c (q,n) 



a c {q) 



[l-exp(-2o c (g)Tn)] , 



(31) 



(32) 



(33) 



(34) 



where T = T\+T 2 is the average duration for one complete cycle with T± = pT and T 2 = (1 — p)T. The definitions 
of a c and D c are same as the uniform duration case. Hence, in terms of the mean durations of the primary processes, 
the above SF is identical to the SF (10) of uniform duration. Thus the scaling exponents corresponding to the mean 
cyclical roughness will be identical to the scaling exponents of the primary process with smaller z, as observed in the 
case of constant duration of primary processes. 

This may be seen in Fig. 1, where two linear processes namely EW and DT (MH universality) growth models 
are used to simulate cyclical growth with random durations distributed uniformly between two values. The scaling 
exponents obtained (a — 0.50 ±0.02 and (3 — 0.25 ±0.03) are EW exponents. However, the mean roughness is slightly 



sinh{2a i (AT i /2)} 



> 1 



higher (Fig. 2) than the roughness of the corresponding uniform duration case. This is expected as — { 2a (AT/2)} 
and hence from eqs. (31) and (32), (£) > cxp(— 2a\T\) and (rj) > exp(— 2a 2 T 2 ). 

If we include nonlinear processes in our simulations the scaling behavior still remain unchanged from the corre- 
sponding uniform duration case. In Fig. 3, the mean roughness is plotted against the number of cycles (for different 
system sizes) using two nonlinear models, namely, the KK model (KPZ universality) and the LD model (MBE uni- 
versality). The durations of those primary processes were taken to be uniformly distributed with identical average. 
The result: a — 0.50 ±0.02 and (i = 0.32 ±0.03 is consistent with the KPZ values same as obtained in the nonrandom 
case. 



The distribution function is given by P(ti) 
distribution function 



B. Poisson Distribution 

exp(-T i )T* i 



, with Tj being the mean. Averaging over this discrete 
(0 = ^M-2at f M -J^ (35) 



t=o 



= exp{Ti(exp(-2oi)-l)}, 



(36) 



where T\ is the mean duration for the first process. Similarly, (77) = cxp(T 2 (cxp(— 2a 2 ) — 1) with T 2 being the 
average duration of the second process. 

As before, to extract the scaling exponents small q, i.e. small di (~ q Zi ), limit is taken. In this limit (£) ~ exp{Ti(l — 
2ai — 1)} = cxp(— 2a\T\), as obtained in the uniform duration case. Hence the scaling exponents corresponding to 
the average roughness will be same as the uniform duration case. 
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C. Exponential Distribution 



The distribution function P{ti) = jr exp(— ti/Ti), where T; is the mean value. This yields 

^ = Trh^ andW = rriby (37) 

with T\ and T 2 being the average duration of the two primary processes. In the asymptotic scaling limit q — > 0, 
(£) ~ 1 — 2aiTi and (77) ~ 1 — 2aiX 2 . Since a, ~ <? Zi , the process with smaller z$ will dominate the SF (19). Hence, 
the scaling exponents are once again identical to those of the non-random case. 

D. Power-law Distribution 

If the long time behavior of a primary process has a fat tail, then duration of that process can be taken to be 
distributed as a power law of time with 

P{t) = /it-("+i) t > 1 

= 0, t<l (38) 

where /j is the only parameter of the distribution. Average duration is easily calculated for the above p. d.f. 

(t) = r p(t)dt = n r t-^ +i ut = |~ (39) 

J 1 J 1 1 — A* 

There are two possibilities: 

(i) For [i > 1: (t) = j/^j- Hence the average and all the higher moments are finite. 

(ii) For fi < 1: (t) — ► 00, the average diverges. 
Below we will consider two possibilities separately. 

1. yt > 1 

In the scaling limit q — > 0, the averages of the random durations £ and 77 are given by (see eq. (23)) 

(0^1- 2aiTi and (77) ~ 1 - 2a 2 T 2 . (40) 

From the expression of the SF (19) we can see that the process with smaller z will dominate the average SF and 
hence the roughness. Hence, the scaling exponents should not change from their corresponding uniform duration case. 

We simulated MH deposition and EW desorption alternately, with their duration distributed with a power law tail 
(p > 1). Fig. 4 shows the scaling of the average roughness with number of cycles for different values of fx > 1. The 
exponent (3 ~ 0.25 for all those cases, which is the (3 of the dominating EW process. The roughness exponent a should 
not change from its uniform duration value and we obtain the EW exponent a (= 0.48 ± 0.05) in our simulations. 

2. /j < 1 

In section III we have already analyzed this case of diverging mean duration. The growth exponent (3 is supposed to 
change continuously depending on the value of the parameter /i. This behavior was observed (Fig. 5) in our simulation 
of the cyclical process combining (1+1 )d EW and MH model. Random durations of both the primary processes obey 
power law distributions with identical \i. The roughness exponent a retains its value (Fig. 6) corresponding to uniform 
application time (i.e. the EW value). 

Presence of nonlinearity in one or both the primary processes did not change the essential scaling behavior from the 
linear case. We simulated nonlinear KK (KPZ) and LD (MBE) models having durations distributed with power law 
tails (both /i > 1 and ji < 1). For ji > 1, no difference was observed in the scaling exponents from the corresponding 
uniform duration case. For fi < 1, however, the growth exponent (3 increased continuously as the parameter [i decreased 
to zero, consistent with the scaling behavior found for the linear primary processes. The roughness exponent a remain 
unchanged in all those cases. 

Therefore, for intermittent applications with a diverging average, the dynamic exponent and the growth exponent 
start to differ from those of the uniform case whereas the roughness exponent remains the same. 
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V. CONCLUSIONS 



We are now in a position to provide the answers to the two central questions we posed at the beginning of this 
paper: 

1. To what extent the cyclical growth behavior describes also randomly alternating processes? - The scaling behavior 
of the perfectly periodic growth processes continues to hold for random application times as long as the latter have 
finite means. Small fluctuations in the application times are thus irrelevant. For linear systems it means that the 
primary process with the smallest dynamic exponent will dominate. For non-linear processes, the approximate RG 
procedure devised for cyclical processes may still be applied to extract the asymptotic scaling behavior. 

2. When cyclical growth fails to describe randomly alternating processes, what are the their new scaling properties? 
If at least the distribution of one of the application time decays as a power law i -1 ^ with [i < 1 (such that its mean 
diverges), the effective dynamic exponent is reduced by a multiplicative factor /i. For linear processes it definitely 
implies that the process with the smallest effective dynamic exponent dominates (the value of fj, = 1 should be assigned 
for processes with a finite mean). We claim that the same continues to hold for non-linear processes as well. This 
claim follows from the fact that a primary process with a diverging mean application time is dominated by a few 
especially long applications and therefore renormalizes independently from the other primary process. Theoretical 
studies to further confirm this behavior will be most welcome. 

As far as experimental investigations are concerned, they should be straightforward for systems with finite means 
(e.g., by methods similar to those used to measure the scaling behavior of cyclical growth). It will be more challenging 
to conduct experimental investigations of surface growth by randomly alternating processes where at least one of them 
has a diverging mean application time. 
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FIG. 1. In IF (roughness) vs In n (number of applications) of the MH/EW process with finite-mean random application times 
[Inset: lnIF s (maximal roughness) vslnL]. 
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FIG. 2. The roughness W vs number of applications n of the MH/EW process with finite-mean randomly alternating 
durations, compared with the uniform (i.e. cyclical) case (log-log plot). 

FIG. 3. In W (roughness) vs Inn (number of applications) of the randomly alternating MBE/KPZ process for different system 
sizes L [Inset: \nW s (maximal roughness) vs InL]. 



FIG. 4. In W vs Inn of the power-law distributed MH/EW randomly alternating process for different values of fi > 1. 



FIG. 5. \nW vs Inn of the power-law distributed MH/EW randomly alternating process for different values of fi < 1. 

FIG. 6. The roughness W vs number of applications n of the power-law (fi = 0.75) distributed MH/EW process for different 
system sizes L (log-log plot) [Inset: In W a (maximal roughness) vs InL]. 
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